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(N 

I Abstract. We state a conjecture relating integer-valued superharmonic functions 

on T? to an Apollonian circle packing of R^. The conjecture is motivated by 
the Abelian sandpile process, which evolves configurations of chips on the integer 
lattice by toppling any vertex with at least 4 chips, distributing one of its chips to 
each of its 4 neighbors. When begun from a large stack of chips, the terminal state 
^ of the sandpile has a curious fractal structure which has remained unexplained. 

Our conjecture implies that the Sandpile PDE recently shown to characterize the 
continuum limit of the sandpile is equivalent to the Apollonian PDE, and we use 
the special geometric structure of the latter to prove that it admits certain fractal 
solutions. Boundary condition evidence from finite sandpiles suggest that these 
solutions exactly correspond to regions of the limiting sandpile, leading to precise 
geometric conjectures on the Abelian sandpile's fractal behavior. 

1- Introduction 

s 

I— —I 1.1. Background. First introduced in 1987 by Bak, Tang and Wiesenfeld [I] as a 

model of self-organized criticality, the Abelian Sandpile is equally fascinating as a model 
" ^ of pattern formation. In its simplest form, the sandpile process evolves a configuration 

77 : — > N of chips by iterating a simple rule: find a lattice point x al? with at least 
four chips and topple it, moving one chip from x to each of its four lattice neighbors. 
OO When the initial configuration has finitely many total chips, the sandpile process 

always finds a stable configuration, where each lattice point has at most three chips. 
00 Dhar [6] observed that the resulting stable configuration does not depend on the top- 

pling order, which is the reason for terming the process "Abelian." When the initial 
configuration consists of a large number of chips at the origin, the final configuration 
has a curious fractal structure [3,8, 15-17] which (after rescaling) is insensitive to the 
number of chips. In 25 years of research (see [ ' '] for a brief survey, and [7,20] for more 
detail) this fractal structure has resisted explanation or even a precise description. 

If s„ : — > N denotes the stabilization of n chips placed at the origin, then the 
rescaled configurations 

Sn{x) := Sn{[n^^^x]) 

(where [x] indicates a closest lattice point to a; G M^) converge to a unique limit Sqo- 
This article presents a partial explanation for the apparent fractal structure of this 
limit. 
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Figure 1 . The boundary of F. The shade of gray at location (a, 5) € 
[0,4] X [0,4] indicates the largest c e [2,3] such that M(a, &, c) S F. 
White and black correspond to c = 2 and c = 3, respectively. 



The convergence s„ — was obtained Pegden-Smart [i^], who used viscosity 
solution theory to identify the continuum limit of the least action principle of Fey- 
Levine-Peres [lU]. We call a 2 x 2 real symmetric matrix A stabilizable if there is a 
function u : — Z such that 

u{x) > ^x^Ax and A'^u{x) < 3, (1.1) 

for all X £ 1? ^ where 

AMx) = Y.^u{y)-u[x)) (1.2) 

is the discrete Laplacian of u on 1? and the sum is taken over the four lattice neigh- 
bors. (We establish a direct correspondence between stabilizable matrices and infinite 
stabilizable sandpile configurations in Section 3.) It turns out that the closure f of the 
set F of stabilizable matrices determines s^o ■ 
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Theorem 1.1 (Existence Of Scaling Limit, [ ]). The resettled configurations Sn con- 
verge weakly-* in L°°(M^) to Srx, — Auoo, where 

Voo ■■= mm{w e C{R^) | u; > -$ and D'^(w + G f }. (1.3) 

Here^{x) -(27r)"Mog| x\ is the fundamental solution of the Laplace equation — 
0, the minimum is taken pointwise, and the differential inclusion is interpreted in the 
sense of viscosity. 

Recall that weak-* convergence simply captures convergence of the local average 
value of Sn- The sum Uqo — v^o + $ has a natural interpretation in terms of the 
sandpile: it is the limit Uao{x) — \\vt\n^^n^^Un{[n^^'^x\)^ where u^ix) is the number 
of times x topples during the formation of s„. As explained below in Section 2.3, 
the function Uao solves the Sandpile PDE, 

D^u e or, (1.4) 

in the open set {u > 0} where it is above the obstacle. 

1.2. Apollonian structure. The key players in the obstacle problem (1.3) are $ and 
r. The former encodes the initial condition (with the particular choice of — (27r)~^ log |a;| 
corresponding to all particles starting at the origin). The set F is a more interesting 
object: it encodes the continuum limit of the sandpile stabilization rule. In Section 3, 
we present an algorithm to determine whether a given matrix with rational entries lies 
in r. The numerical evidence provided by this algorithm indicates that F is a union of 
downward cones based at points of a certain set V. The elements of V, which we call 
peaks, are visible as the locally darkest points in Figure 1. 

We can state an exact conjecture characterizing the the shape of F in terms of Apol- 
lonian configurations of circles. Three pairwise externally tangent circles Ci,C2,C3 
determine an Apollonian circle packing, as the smallest set of circles containing them 
that is closed under the operation of adding, for each pairwise tangent triple of circles, 
the two circles which are tangent to each circle in the triple. They also determine a 
downward Apollonian packing, closed under adding, for each pairwise-tangent triple, 
only the smaller of the two tangent circles. Lines are allowed as circles, and the Apollo- 
nian band circle packing is the packing C determined by the lines {x = 0} and {x = 2} 
and the circle {{x — 1)^ + y'^ — !}• Its circles are all contained in the strip [0, 2] x M. 

We put the proper circles in (i.e., the circles that arc not lines) in bijective 
correspondence with real symmetric 2x2 matrices of trace > 2, in the following way. 
To a proper circle C = {(a; — a)^ + (y — 6)^ = r^} in we associate the matrix 

m(C) :=M(a,&,r-f2) 

where 

' +a b 

b c- 



M{a,b,c) -.^ i 



(1.5) 



We write 5*2 for the set of symmetric 2x2 matrices with real entries, and, for A,BgS2 
we write B < A ii A — B is nonnegative definite. For a set C 52, we define 

:= {B e 82] B < A for some A e V}, 

the order ideal generated by V in the matrix order. 

Now let C = Ufeez(^ + (2^' *-*)) be the extension of the Apollonian band packing to 
all of E^ by translation. Let 

V = {m{C) I C e C}. 
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Figure 2. An Apollonian triangulation is a union of Apollonian 
triangles meeting at right angles, whose intersection structure matches 
the tangency structure of their corresponding circles. The solution u 
of Theorem 1.3 has constant Laplacian on each Apollonian triangle, 
as indicated by the shading (darker regions are where Au is larger). 

Conjecture 1.2 (Circle Packing Conjecture). F = V^. 

Note that replacing the condition A^v < 3 with A^v < in (1.1) results only in a 
translation of the set F. Indeed, on can simply subtract the function v{x) — |xi(a;i + l) 
from u and subtract Af(3,0, 3) from A. Thus this conjecture can be interpreted as a 
statement on the geometry of the set of integer-valued superharmonic functions of 
quadratic growth on the lattice Z^. The set of 4-tuples of pairwise tangent circles 
in an Apollonian circle packing has a transitive action by a discrete subgroup of the 
Lorenz group 0(3, 1), and group theory has enabled recent breakthroughs in proving 
deep arithmetic properties of Apollonian packings [2, 11, 19]. Conjecture 1.2 has so far 
resisted any attempts to bring these tools to bear. 

1.3. The Apollonian PDE. Conjecture 1.2 asserts that the Sandpile PDF (1.4) is 
the same as the Apollonian PDE 

D^u e dV^. (1.6) 

Our main result. Theorem 1.3 below, constructs a family of piecewise quadratic solu- 
tions to the Apollonian PDF. The supports of these solutions are the closures of certain 
fractal subsets of which we call Apollonian triangulations, giving an explanation 
for the fractal limit Soo- 

Of course, every matrix A £ S2 with tr(j4) > 2 is now associated to a unique proper 
circle C = c{A) = m~^{A) in K^. We say two matrices are {externally) tangent pre- 
cisely if their corresponding circles are (externally) tangent. Given pairwise externally 
tangent matrices Ai,A2, A3, denote by ^(^1, A2, A3) (resp. {Ai, A2, A^)) the set of 
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Figure 3. Left: The sandpile s„ for n = 4 • 10^. Sites with 0, 1, 2, 
and 3 chips are represented by four different shades of gray. Right: 
A zoomed view of the boxed region, one of many that we beheve 
converges to an ApoUonian triangulation in the n ^ oo hmit. 



matrices corresponding to the Apollonian circle packing (resp. downward Apollonian 
packing) determined by the circles corresponding to ^1,^2,^3. 

Theorem 1.3 (Piecewise Quadratic Solutions). For any pairwise externally tangent 
matrices ^1,^2,^3 € S2, there is a nonempty convex set Z C M."^ and a function 
u £ C^'^(Z) satisfying 

in the sense of viscosity. Moreover, Z decomposes into disjoint open sets (whose clo- 
sures cover Z) on each of which u is quadratic with Hessian in A~ {Ai^ A2, A'i) . 

This theorem is illustrated in Figure 2. We call the configuration of pieces where D^u 
is constant an Apollonian triangulation. Our geometric characterization of Apollonian 
triangulations begins with the definition of Apollonian curves and Apollonian triangles 
in Section 5. We will see that three vertices in general position determine a unique 
Apollonian triangle with those vertices, via a purely geometric construction based on 
medians of triangles. We will also show that any Apollonian triangle occupies exactly 
4/7 of the area of the Euclidean triangle with the same vertices. 

An Apollonian triangulation, which we precisely define in Section 6, is a union of 
Apollonian triangles corresponding to circles in an Apollonian circle packing, where 
pairs of Apollonian triangles corresponding to pairs of intersecting circles meet at 
right angles. The existence of Apollonian triangulations is itself nontrivial and is the 
subject of Theorem 7.1; analagous discrete structures were constructed by Paoletti 
in his thesis [17]. Looking at the Apollonian fractal in Figure 2 and recalling the 
51/2 (2) symmetries of Apollonian circle packings, it is natural to wonder whether 
nice symmetries may relate distinct Apollonian triangulations as well. But we will 
see in Section 6 that Apollonian triangles are equivalent under affine transformations, 
precluding the possiblity of conformal equivalence for Apollonian triangulations. 
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If Ci, C2, C3 are pairwise tangent circles in the band circle packing, then letting Ai — 
m{Ci) for 1 = 1,2, 3, we have A~ {Ai, A2, A^) C V, so the function u in Theorem 1.3 will 
be a viscosity solution to the Sandpile PDE assuming Conjecture 1.2. The uniqueness 
machinery for viscosity solutions gives the following corollary to Theorem 1.3, which 
encapsulates its relevance to the Abelian Sandpile. 

Corollary 1.4. Suppose Ui,U2,Uj, C are connected open sets bounding a convex 
region Z such that Ui Ci Uj = {xk} for {«, j, k} = {1, 2, 3}, where the triangle Axia;2a;3 
is acute. If u^o is quadratic on each of Ui,U2,U3 with pairwise tangent Hessians 
Ai, A2, A3 G V, respectively, then Conjecture 1.2 implies that u^o is piecewise qua- 
dratic in R and the domains of the quadratic pieces form the Apollonian triangulation 
determined by the vertices Xi,X2,X3. 

Note that Soo = Auoo implies Soo is piecewise-constant in the Apollonian triangulation. 

Let us briefly remark on the consequences of this corollary for our understanding 
of the limit sandpile. As observed in [S, 111] and visible in Figure 3, the sandpile 
s„ for large n features many clearly visible patches, each with its own characteristic 
periodic pattern of sand (sometimes punctuated by one-dimensional 'defects' which 
are not relevant to the weak-* limit of the sandpile). Empirically, we observe that 
triples of touching regions of these kinds are always regions where the observed finite w„ 
correspond (away from the one-dimensional defects) exactly to minimal representatives 
in the sense of (1.1) of quadratic forms 

^x*Ax + bx 

where the A's for each region are always as required by Corollary 1.4. Thus we are 
confident from the numerical evidence that the conditions required for Corollary 1.4 and 
thus Apollonian triangulations occur — indeed, are nearly ubiquitous — in Sao- Going 
beyond Corollary 1.4's dependence on local boundary knowledge would seem to require 
an understanding the global geometry of Soo, which remains a considerable challenge. 

1.4. Overview. The rest of the paper proceeds as follows. In Section 2, we review 
some background material on the abelian sandpile and viscosity solutions. Section 
3 presents our algorithm for verifying Conjecture 1.2 to arbitrary precision. After 
reviewing some basic geometry of Apollonian circle packings in Section 4, we define 
and study Apollonian curves, Apollonian triangles, and Apollonian triangulations in 
Sections 5 and 6. The proofs of Theorem 1.3 and Corollary 1.4 come in Section 7 where 
we construct piecewise-quadratic solutions to the Apollonian PDE. Finally, in Section 
8 we discuss new problems suggested by our results. 

2. Preliminaries 

The preliminaries here are largely section-specific, with Section 2.1 being necessary 
for Section 3 and Sections 2.2 and 2.3 being necessary for Section 7. 

2.1. The Abelian sandpile. Given a configuration 77 : ^ Z of chips on the integer 
lattice, we define a toppling sequence as a finite or infinite sequence xi, X2, X3, . . . of 
vertices to be toppled in the sequence order, such that any vertex topples only finitely 
many times (thus giving a well-defined terminal configuration). A sequence is legal if 
it only topples vertices with at least 4 chips, and stabilizing if there are at most 3 chips 
at every vertex in the terminal configuration. We say that 77 is stabilizable if there 
exists a legal stabilizing toppling sequence. 
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The theory of the Abehan Sandpile begins with the foUowing standard fact: 

Proposition 2.1. Any x ^1? topples at most as many times in any legal sequence 
as it does in any stabilizing sequence. □ 

Proposition 2.1 imphes that to any stabihzable initial configuration rj, we can associate 
an odometer function u : — > N which counts the number of times each vertex topples 
in any legal stabilizing sequence of topplings. The terminal configuration of any such 
sequence of topplings is then given hy ij + A^v. Since v and so A^v are independent of 
the particular legal stabilizing sequence, this shows that the sandpile process is indeed 
"Abelian" : if we start with some stabilizable configuration rj > 0, and topple vertices 
with at least 4 chips until we cannot do so any more, then the final configuration 
T] + A^v is determined by rj. 

The discrete Laplacian is monotone, in the sense that A^u{x) is decreasing in u{x) 
and increasing in u{y) for any neighbor ?/ ~ x of a; in Z^. An obvious consequence of 
monotonicity is that taking a pointwise minimum of two functions cannot increase the 
Laplacian at a point: 

Proposition 2.2. If u,v : Z"^ ^ Z, w := min{u, v}, and w{x) — u{x), then A^w{x) < 
A^u{x). □ 

In particular, given any functions u, v satisfying rj + A^{u) < 3 and rj + A^{v) < 3, 
their pointwise minimum satisfies the same constraint. The proof of Theorem 1.1 in 
[li^] begins from the Least Action Principle formulated in [Id], which states that the 
odometer of an initial configuration rj is the pointwise minimum of all such functions. 

Proposition 2.3 (Least Action Principle). If rj : 1? ^ 1^ and w : 2? ^ 1^ satisfy 
rj + A^w < 3, then rj is stabilizable, and its odometer v satisfies v < w. 

Note that the Least Action Principle can be deduced from Proposition 2.1 by as- 
sociating a stabilizing sequence to w. By considering the function u = v — 1 for any 
odometer function v, the Least Action Principle implies the following proposition: 

Proposition 2.4. If rj : "Z"^ ^ "Z is a stabilizable configuration, then its odometer v 
satisfies v{x) = for some x . 

Finally, we note that these propositions generalize in a natural way from Z? to 
arbitrary graphs; in our case, it is sufhcient to note that they hold as well on the torus 

T„:=ZV«Z2 forneZ+. 

2.2. Some matrix geometry. All matrices considered in this paper are 2x2 real 
symmetric matrices and we parametrize the space 5*2 of such matrices via M : M'^ — )■ S'2 
defined in (1.5). We use the usual matrix ordering: A < _B if and only if _B — A is 
nonnegative definite. 

Of particular importance to us is the downward cone 

A^ ■= {B S2 ■■ B < A}. 

Recall that if i? S dA^ , then A — B — v ® v — for some column vector v. That is, 
the boundary dA^ consists of all downward rank-1 perturbations of A. 

Our choice of parametrization M was chosen to make A^ a cone in the usual sense. 
Observe that 

M(a,6,c)>0 if and only if c> (c? ^b^fl"^ . 

Moreover: 
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Observation 2.5. We have 

v(S)v = M{ui,U2,{ul+uiy^^) (2.1) 
if and only if v"^ — u, where v"^ denotes the complex square of v. □ 
Thus if S e dA^, then 

A-B = {-p{B) - ^ _ -p[A)Y/\ (2.2) 

where 

-p{M{a,h,c)) := (a, 5), 
and w^/^ denotes the complex square root of a vector w e = C. 
Denoting by / the 2x2 identity matrix, we write 

A^ = A- 2(tr(A) - 2)1 

for the reflection of A across the trace-2 plane; and 

^0^^ + ^- 



2 

for the projection of A on the trace-2 plane. Since the line {A + t{v (8) w) | ^ G M} is 
tangent to the downward cone A^ for every nonzero vector v and matrix A, we see that 
matrices Ai,A2, both with trace greater than 2, are externally tangent if and only if 
Ai — A2 has rank 1 and internally tangent if and only if Ai — A2 has rank 1. This 
gives the following Observation: 

Observation 2.6. Suppose the matrices Ai,Aj,Ak are mutually externally tangent 
and have traces > 2. Then there are at most two matrices B whose difference As — B 
is rank 1 for each s = i,j,k: B — A^ is a solution for any matrix A^ externally 
tangent to Ai, Aj, Ak, and B = A„i is a solution for any A„i internally tangent to 
Ai , Aj , Ak ■ O 

Note that the case of fewer than two solutions occurs when the triple of trace-2 circles 
of the down-set cones of the Ai are tangent to a common line, leaving only one proper 
circle tangent to the triple. 

2.3. Viscosity Solutions. We would like to interpret the Sandpile PDE D'^u G dT in 
the classical sense, but the nonlinear structure of dT makes this impractical. Instead, 
we must adopt a suitable notion of weak solution, which for us is the viscosity solution. 
The theory of viscosity solutions is quite rich and we refer the interested reader to [4,5] 
for an introduction. Here we simply give the basic definitions. We remark that these 
definitions and results make sense for any subset F C 5*2 that is downward closed and 
whose boundary has bounded trace (see Facts 3.2, 3.5, and 3.6 below). 

If C is an open set and u S C(0), we say that u satisfies the differential 
inclusion 

D'^u e f in rj, (2.3) 
if D^ip{x) G F whenever ip € C°°(ri) touches u from below a,t x £ H.. Letting F^ denote 
the closure of the complement of F, we say that u satisfies 

D^u e F^ in (2.4) 

if D'^ip{x) e F^ whenever t/j e C°°{il) touches u from above a,t x E fl. Finally, we say 
that u satisfies 

D^u e dV in O, 

if it satisfies both (2.3) and (2.4). 
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The standard machinery for viscosity solutions gives existence, uniqueness, and 
stabiUty of solutions of the Sandpile PDE. For example, the minimum in (1.3) is 
indeed attained by some v G C(]R^) and u + $ solves the Sandpile PDE in {v > — $}. 
Moreover, we have a comparison principle: 

Proposition 2.7. //il C is open and bounded and u,v G C(f2) satisfy 

D^u e f and D'^v e m n, 

then supj2(w — u) — supgQ(t; — u). □ 

Recall that C^'^{U) is the class of differentiable functions on U with Lipshitz deriva- 
tives. In Section 7, we construct piecewise quadratic C^'^ functions which solve the 
Apollonian PDE on each piece. The following standard fact guarantees that the func- 
tions we construct are, in fact, viscosity solutions of the Apollonian PDE on the whole 
domain (including at the interfaces of the pieces). 

Proposition 2.8. // J7 C is open, u G C^'^{U), and for Lebesgue almost every 
xGU 

D^u{x) exists and D^u{x) e dT, 
then D^u G dT holds in the viscosity sense. □ 

Since we are unable to find a published proof, we include one here. 

Proof. Suppose G C°°{U) touches u from below at xq G U. We must show D'^(f{xo) G 
r. By approximation, we may assume that is a quadratic polynomial. Fix a small 
£ > 0. Let A be the set oi y G U for which there exists p G and g G M such that 

(py{x) := (p{x) - ^e\x\'^ +p- x + q, 

touches u from below a y. Since u G C^'^, p{y) is unique and that map p : A ^ M? 
is Lipschitz. Since e > and U is open, the image p{A) contains a small ball -65(0). 
Thus we have 

0< 15^(0)1 < \p{A)\ < Lip{p)\A\. 
In particular, A has positive Lebesgue measure and we may select a point y G A such 
that D'^u{y) exists and D'^u{y) G f. Since (ph touches u from below at y, we have 
D'^fyiu) < D'^u{y) and thus D'^ipy(y) = D'^tp{y) - el = D'^ip(xo) - el gT. Sending 
e 0, we obtain D^(p{xo) G f . □ 

3. Algorithm to decide membership in F 

The definition of F does not give a method for verifying membership in the set. In 
this section, we will show that matrices in F correspond to certain infinite stabilizable 
sandpiles Z^. li A G T has rational entries, then its associated sandpile is periodic, 
which yields a method for checking membership in F for any rational matrix, and 
allows us to algorithmically determine the height of the boundary of F at any point 
with arbitrary precision. If g : — >■ K, write \q] for the function — > Z obtained 
by rounding each value of q up to the nearest integer. The principal lemma is the 
following. 

Lemma 3.1. A G T if and only if the configuration \qA] is stabilizable, where 

qA{x) := ^x^Ax 
is the quadratic form associated to A. 
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Proof. If u satisfies (f .1), tlicn tlie Least Action Principle applied to w = u — Iqa] 
shows that 77 = Iqa] is stabilizable. On the other hand, if r/ = Iqa] is stabilizable 
with odometer v, then u — v + Iqa] satisfies (f .1). □ 

Since A < B implies x^Ax < x^Bx for all x e Z^, the definition of F implies that F 
is downward closed in the matrix order: 

Fact 3.2. IfA<B and B eT, then AeT. □ 

It follows that the boundary of F is Lipschitz, and in particular, continuous; thus to 
determine the structure of F, it suffices to characterize the rational matrices in F. We 
will say that a function s on is n-periodic if s{x + y) ~ s{x) for all y e nZ^. 

Lemma 3.3. // A has entries in for a positive integer n, then Iqa] is 'in- 
periodic. 

Proof. If 1/ e 2nl? then Ay e 21?, so 

qA{x + y) - qA{x) = (a;* + \^y'-)Ay G Z. 
Hence \qA\ — qa is 2n-periodic. Writing 

Al \qA^ = A^l qA^ - q a) - A\a 
and noting that A^qA is constant, we conclude that A^ \qA^ is 2n-periodic. □ 

Thus the following lemma will allow us to make the crucial connection between rational 
matrices in F stabilizable sandpiles on finite graphs. It can be proved by appealing to 
[9, Theorem 2.8] on infinite toppling procedures, but we give a self-contained proof. 

Lemma 3.4. An n-periodic configuration rj : 7? ^ Z is stabilizable if and only if it is 
stabilizable on the torus r„ — 1? jril? . 

Proof. Supposing 77 is stabilizable on the torus T„ with odometer w, and extending v 
to an n-periodic function v on in the natural way, we have that 77 -f A^v < 3. Thus 
?7 is stabilizable on Z^ by the Least Action Principle. 

Conversely, if r] is stabilizable on Z^ , then there is a function w : Z'^ — )■ N such that 
r] -\- A^w < 3. Proposition 2.2 implies that 

w{x) :— mm{w{x + y) '■ y & nl?}, 

also satisfies 77 -I- A^w < 3. Since w is 7i-periodic, we also have 77 -I- A^p w < 3 and thus 
77 is stabilizable on the torus T„. □ 

The preceding lemmas give us a simple prescription for checking whether a rational 
matrix A is in F: compute s — A^ \qA^ on the appropriate torus, and check if this 
is a stabilizable configuration. To check that s is stabilizable on the torus, we simply 
topple vertices with > 4 chips until either reaching a stable configuration, or until 
every vertex has toppled at least once, in which case Proposition 2.4 implies that s is 
not stabilizable. 

We thus can determine the boundary of F to arbitrary precision algorithmically. 
For (a, 6) e let us define 

co(a, b) = sup{c I M (a, 6, c) e F}. 

By Fact 3.2, we have M{a,b,c) G f if and only if c < co{a,b). Hence the boundary 
dT is completely determined by the Lipschitz function Co(a, b). In Figure 1, the shade 
of the pixel at (a, 5) corresponds to a value c that is provably within of Co(a,6). 
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(Note that because of the Lipshitz condition, for any e > there is a finite procedure 
to either prove that T is within e of its conjectured shape, or disprove the conjecture.) 

The above resuhs ahow for some very Umited partial confirmations of Conjecture 1.2. 
In particular, it is easy to deduce the following two facts: 

Fact 3.5. If A is rational and tr{A) < 2, then A eT. □ 
Fact 3.6. If A is rational and tr{A) > 3, then A^T. □ 

In both cases, the relevant observation is that for rational A, tr{A) is exactly the 
average density of the corresponding configuration 77 = Iqa] on the appropriate 
torus. This is all that is necessary for Fact 3.6. For Fact 3.5, the additional ob- 
servation needed (due to Rossin [21]) is that on any finite connected graph, a chip 
configuration with fewer chips than there are edges in the graph will necessarily stabi- 
lize: for unstabilizable configurations, a legal sequence toppling every vertex at least 
once gives an injection from the edges of the graph to the chips, mapping each edge to 
the last chip to travel across it. 

Facts 3.5 and 3.6 along with continuity imply that 2 < co(a, 6) < 3 for all (a, b) E M^. 
With some additional work, the above results can be used to show that co{a,b) = 2 
for all a G 2Z and 6 S M, confirming Conjecture 1.2 along the vertical lines x — a for 
a e 2Z. Finally, let us remark that cq has the translation symmetries 

Co(a + 2, 6) = co(a, b) = co(a, b + 2). 

This follows easily from the observation that ^x{x + 1) — ^y{y + 1) and xy are both 
integer- valued discrete harmonic functions on Z^. 

4. Apollonian circle packings 

For any three tangent circles Ci,C2,C3, we consider the corresponding triple of 
tangent closed discs Di, D2, D3 with disjoint interiors. We allow lines as circles, and 
allow the closure of any connected component of the complement of a circle as a closed 
disc. Thus we allow internal tangencies, in which case one of the closed discs is actually 
the unbounded complement of an open bounded disc. Note that to consider Ci, C2, C3 
pairwise tangent we must require that three pairwise intersection points of the Ci are 
actually distinct, or else the corresponding configuration of the Di is not possible. In 
particular, there can be at most two lines among the Ci, which are considered to be 
tangent at infinity whenever they are parallel. 

The three tangent closed discs Di,D2, D3 divide the plane into exactly two regions; 
thus any pairwise triple of circles has two Soddy circles, tangent to each circle in the 
triple. If all tangencies are external and at most one of Ci , C2 , C3 is a line, then 
exactly one of the two regions bordered by the Di is bounded, and the Soddy circle in 
the bounded region is called the successor of the triple. 

An Apollonian circle packing, as defined in the introduction, is a minimal set of cir- 
cles containing some triple of pairwise-tangent circles and closed under adding all Soddy 
circles of pairwise-tangent triples. Similarly, a downward Apollonian circle packing is 
a minimal set of circles containing some triple of pairwise externally tangent circles 
and closed under adding all successors of pairwise-tangent triples. 

For us, the crucial example of an Apollonian packing is the Apollonian band packing. 
This is the packing which appears in Conjecture 1.2. A famous subset is the Ford 
circles, the set of circles Cp/g with center ( y , ^ ) and radius ^ , where p/q is a rational 
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number in lowest terms. A simple description of the other circles is, as far as we know, 
unavailable. Thus our conjecture provides a new perspective on these circles. 



An important observation regarding Apollonian circle packings is that a triple of 
pairwise externally tangent circles is determined by its intersection points with its 
successor: 

Proposition 4.1. Given a circle C and points 2/1,2/2,2/3 S C, there is at exactly one 
choice of pairwise externally tangent circles Ci,C2,C3 which are externally tangent to 
C at the points yi,y2,y3- O 

Proposition 4.1, together with its counterpart for the case allowing an internal tan- 
gency, allows the deduction of the following fundamental property of Apollonian circle 
packings. 

Proposition 4.2. Let C be an Apollonian circle packing. A set C of circles is an 
Apollonian circle packing if and only if C = fJ-{C) for some Mohius transformation 
fi. □ 

The use of Mobius transformations allows us to deduce a geometric rule based on 
medians of triangles concerning successor circles in Apollonian packings: 

Lemma 4.3. Suppose that circles C, Ci,C2 are pairwise tangent, with Soddy circles 
Co and C3, and let zf = pi — c, viewed as a complex number, where c is the center of 
C and Pi is the intersection point of C and Ci for each i. If Li is a line parallel to the 
vector Zi which passes through if i ~ 1, 2, 3 and does not pass through if i = 0, then 
Lt, is a median line of the triangle formed by the lines Lq, Li, L2. 

Proof. Without loss of generality, we assume that C is a unit circle centered at the 
origin, and that Zg = —1. The Mobius transformation 

_ zi + Z1Z2 - z{zi - Z2) 

1 + Z2 + z{Zi - Z2) 

sends to z^, 1 to z|, and cx) to — 1 = Zq. Thus, for the pairwise tangent generalized 
circles C = {y = 0},C^ = {2/ = 1}.C[ = {x^Hy-\Y = W.G'2 = {[x-\f^{y^\f = 
j},C3 = {{x — l)'^ = jj^} (these are some of the "Ford circles"), we have that /i maps 
the intersection point of C',C'i to the intersection point of C,Ci for i ^ 0, 1,2, thus 
it must map the intersection point of C", Cg to the intersection point of C, C3, giving 
A'zi,22(|) = ^3- Thus it suffices to show that for 

ft \ n /o\ 2:1 + Z2 + 2ziZ2 

/(zi,Z2) := ^^21,22(1/2) 



we have that 



Zl + Z2 



Re(zi)Im(z2) + Re(z2)Im(zi) . 

— t 



/ Re(zi)Im(z2) + Re(z2)Im(zi) \ 



1 



V 2Re(zi)Re(z2) / 



as the right-hand side is the square of the unit vector whose tangent is the average of the 
tangents of Zi and Z2; this is the correct slope of our median line since Zg = — 1 implies 
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Figure 4. The circle arrangement from Proposition 4.5. 

that Lq is vertical. We will check (4.1) by writing zi — cos a+i sin a, Z2 — cos /3+i sin /3 
to rewrite f{zf, z|) as 



(cos a + isina)^ + (cos/3 + zsin/3)^ + 2(cosa + isina)^(cos/3 + isin/?)^ 

(4.2) 



(cos a + i sin a)^ + (cos /3 + i sin /3)2 + 2 
(cos(a + /3) + isin(a + /3))(cos(a — /3) + cos(a + /3) + i sin(a + /?)) 



cos(q! — /3)(cos(a + /3) + i sin(a + /?)) + 1 

where we have used the identity 

(cosx + isinx)^ + (cosy + isiny)^ = 2cos(a: — y){cos{x + y) + i sin(a; + y)), 

which can be seen easily geometrically. Dividing the top and bottom of the right side 
of (4.2) by cos(a + (3) + ism{a + (3) gives 

2 2\ cos(a — /?) + cos(a + /3) + i sin(a + /?) 
^^^^'^'^> ~ cos(a -13) + cos(a + P) - i sin(a + (3) ' 

Thus to complete the proof, note that the right-hand side of (4.1) can be can simplified 

\ 2 

-| I cos a sin /3+cos /3 sin ck ■ \ 2 

^ 2 cos a cos /3 *J (cos(a + /?) + cos(q; — (3) + i sin(a + /?)) 
I , ( cos g sin +cos sin a \ ^ (cos(a + /?) + cos(a — (3)y + sin^ (a + 13) 

\ 2 cos a cos ^ y 

cos(a + f3) + cos(a — /3) + i sin(a + f3) 
cos(a + /3) + cos(a — 13) — i sin(a + /3) 

by multiplying the top and bottom by (2 cos a cos /3)^ and using the Euler identity 
consequences 

2 cos a cos (3 = cos(a + /3) — cos(q; — /3) 
cosasin/3 + cos/^sina = sin(Q: + /3). □ 

Remark 4.4. By Proposition 4.2, a set of three points {xi,X2,X3} on a circle C 
uniquely determine three other points {j/i, 2/2,2/3} on C, as the points of intersection of 
C with successor circles of triples {C,Ci, Cj}, where Ci, C2, C3 are the unique triple of 
circles which are pairwise externally tangent and externally tangent to C at the points 
Xi. Since the median triangle of the median triangle of a triangle T is homothetic to 
T, Lemma 4.3 implies that this operation is an involution: the points determined by 
{2/1)2/2)2/3} in this way is precisely the set {xi, a;2, Xa}. 
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We close this section with a coUection of simple geometric constraints on arrange- 
ments of externally tangent circles (Figure 4) , whose proofs are rather straightforward: 

Proposition 4.5. Let Cq,C'i,C2 be pairwise externally tangent proper circles with 
successor C3, and let C4 and C5 be the successors of Cq, Ci, C3 and Cq, C2, C3, respec- 
tively. Letting Ci denote the center of the circle Ci, we have the following geometric 
bounds: 

(1) c,C3Cj <TTfor{i,j} C {0,1,2}. 

(2) Zc4CoC3,Zc5CoC3 < f. 

(3) ZC4C0C3 > IZC5C0C3 (and vice versa). 

(4) ZC4C3C5 > 2 • arctan(3/4). □ 



We build up to Apollonian triangles and triangulations by defining the Apollonian 
curve associated to an ordered triple of circles. This will allow us to define the Apol- 
lonian triangle associated to a quadruple of circles, and finally the Apollonian trian- 
gulation associated to a downward packing of circles. We will define these objects 
implicitly, and then show that they exist and are unique up to translation and ho- 
mothety (i.e., any two Apollonian curves 7,7' associated to the same triple satisfy 
7' = 07 + b for some a G M and 6 S M^). In Section 6, we give a recursive description 
of the Apollonian curves which characterizes these objects without reference to circle 
packings. 

Fix a circle Cq with center cg and let C and C be tangent circles tangent to Cq 
at X and x' , and have centers c and c', respectively. We define s(C, C") to be the 
successor of the triple (Cq, C, C") and a(C) to be the angle of the vector v{C) := c—cq 
with the positive cc-axis. Let v^^^{C) to be a complex square root of w(C), and let 
^1/2 (c) ^ Mi,i/2((7) be the real hue it spans. (We wiU actually only use ^^/^(C), so 
the choice of square root is immaterial.) Note that all of these functions depend on the 
circle Cq; we will specify which circle the functions are defined with respect to when 
it is not clear from context. 

Now fix circles C\ and C2 such that Cq, Ci, C2 are pairwise externally tangent. Let 
C denote the smallest set of circles such that Ci , C2 € C and for all tangent C,C' ^ C 
we have s(C, C) e C. Note that all circles in C are tangent to Cq. 

Definition 5.1. A (continuous) curve 7 : [q;(Ci), 0(6*2)] is an Apollonian curve 

associated to the triple {Cq, Ci, C2) if for all tangent circles C,C' G C, 



We call 7(a(s(Ci, C2))) the splitting point of 7. The following Observation implies, 
in particular, that the splitting point divides 7 into two smaller Apollonian curves. 

Observation 5.2. For any two tangent circles C,C' G C, the restriction 7|[q(c),q(c)] 



To prove the existence and uniqueness of Apollonian curves, we will need the fol- 
lowing observation, which is easy to verify from the fact that no circle lying inside the 
region bounded by Cq, Ci, C2 and tangent to Cq has interior disjoint from the family 



5. Apollonian triangles and triangulations 



7(a(C))-7(a(C'))e£i/2(s(C,C')). 



is also an Apollonian curve. 



□ 



C: 



Observation 5.3. a(C) is dense in the interval [a(Ci), a(C2)] . 



□ 
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We can now prove the existence and uniqueness of Apollonian curves. 

Theorem 5.4. For any pairwise tangent ordered triple of circles (Cq, Ci, C2), there is 
an associated Apollonian curve 7, which is unique up to translation and scaling. 

Proof. The choice of the points 7(q;(Ci)) and 7(q;(C2)) is determined uniquely up to 
translation and scaling by the constraint that 7(a(Ci)) —j{a{C2)) is a real multiple of 
v^^'^{s{Ci,C2))- This choice then determines the image ^{a{C)) for all circles C € C 
recursively: for any tangent circles C^,C^ G C with s{C^ , C^) the constraints 

j{a{C'))~jia{C^))ef/'{s{C\C')) 

determine 7(a(C'^)) uniquely given 7(a(C^)) and j{a{C'^)). To show that there is a 
unique and well-defined curve 7, by Observation 5.3 it is enough to show that 7 is a 
continuous function on the set a{C). For this it suffices to find an absolute constant 
/3 < 1 for which 

\j{a{C')) - j{ais{C\C')))\ < f3 \i{a{C')) - l{a{C^))\ (5.1) 

for tangent circles C^,C^ G C, as this implies, for example, that by taking successors 
k times, we can find a circle C" G C such that all points in 7([a(C^), q;(C")]) lie 
within ^'^ \^{a{C^)) - 7(a(C^))| of 7(q;(C^)). We get the absolute constant /3 from 
an application of the law of sines to the triangle with vertices pi = ^{a{C^)),p2 = 
j{a{C^)),P3 = -f{a{s{C^, C^))): part 3 of Proposition 4.5 implies that 9 := Zp3pip2 > 
^Zp3P2Pi', the Law of Sines then implies that line (5.1) holds with [3 = , which 

is < I always since part 2 of Proposition 4.5 implies that < f . □ 

Theorem 5.5. The image of an Apollonian curve 7 corresponding to (Co, Ci, C2) has 
a unique tangent line at each point 7(0;). This line is at angle a/2 to the positive 
X-axis. In particular, j is a convex curve. 

Proof. Observation 5.3 and Definition 5.1 give that for any C G C, there is a unique line 
tangent to the image of 7 at 7(a(C)), which is at angle a{C)/2 to the x-axis. Together 
with another application of Observation 5.3 and the fact that ^ is a continuous function 
of a, this gives that the image 7 has a unique tangent line at angle ^ to the x-axis at 
any point 7(0;). □ 

Definition 5.6. The Apollonian triangle corresponding to an unordered triple of ex- 
ternally tangent circles Ci, C2, C3 and circle Cq externally tangent to each of them 
is defined as the bounded region (unique up to translation and scaling) enclosed 
by the images of the Apollonian curves 712, 723, 731 corresponding to the triples 
(Co,Ci,C2), {Co, €2,03), (Co,C3,Ci) such that -fij{a{Cj)) = -fjk{a{Cj)) for each 
{z,i,fc} = {1,2,3}. 

Note that Theorem 5.4 implies that each triple {Ci,C2,C3} of pairwise tangent 
circles corresponds to an Apollonian triangle T which is unique up to translation 
and scaling. Theorem 5.5 implies that the curves 712,723,731 do not intersect ex- 
cept at their endpoints, and that T is strictly contained in the triangle with vertices 
7i2(C2),723(C3),73i(Ci). Another consequence of Theorem 5.5 is that any two sides 
of an Apollonian triangle have the same tangent line at their common vertex. Thus, 
the interior angles of an Apollonian triangle are 0. 
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An Apollonian triangle is proper if C4 is smaller than each of Ci, C2, C3, i.e., if C4 
is the successor of Ci,C2,C3, and all Apollonian triangles appearing in our solutions 
to the Apollonian PDE will be proper. 

We also define a degenerate version of an Apollonian triangle: 

Definition 5.7. The degenerate Apollonian triangle corresponding to the pairwise 
tangent circles (Ci, C2, C3) is the compact region (unique up to translation and scaling) 
enclosed by the image of the Apollonian curve 7 corresponding to (Ci, 6*2,(73), and 
the tangent lines to 7 at its endpoints 7(a(C2)) and 7(a(C3)). 

Proper Apollonian triangles (and their degenerate versions) are the building blocks 
of Apollonian triangulations, the fractals that support piecewise-quadratic solutions 
to the Apollonian PDE. Recall that A^{Ci, C2, C3) denotes the smallest set of circles 
containing the circles Ci, C2, C3 and closed under adding successors of pairwise tangent 
triples. To each circle C € ^~(Ci,C2,C3) \ {Ci,C2,C3} we associate an Apollonian 
triangle Tc corresponding to the unique triple {C^,C^,C^} in {Ci, 02,03) whose 
successor is O. 

Definition 5.8. The Apollonian triangulation associated to a triple {Ci, (72,6*3} of 
externally tangent circles is a union of (proper) Apollonian triangles Tc corresponding 
to each circle O G A~{Oi, O2, O3) \ {Oi, O2, C3}, together with degenerate Apollonian 
triangles Tc for each O = Ci, (72,6*3, such that disjoint circles correspond to disjoint 
Apollonian triangles, and such that for tangent circles 0,0' in A^{Oi,02,03) where 
^{0') < r{0), we have that 7c" and Tc intersect at a vertex of Tc, and that their 
boundary curves meet at right angles. 

Figure 2 shows an Apollonian triangulation, excluding the three degenerate Apollonian 
triangles on the outside. 

Remark 5.9. By Theorem 5.5 and the fact that centers of tangent circles are separated 
by an angle tt about their tangency point, the right angle requirement is equivalent 
to requiring that the intersection of Tc and Tc occurs at the point 7(a((7')) on an 
Apollonian boundary curve 7 of Tc ■ 

6. Geometry of Apollonian curves 

In this section, we will give a circle-free geometric description of Apollonian curves. 
This will allow us to easily deduce geometric bounds necessary for our construction of 
piecewise-quadratic solutions to work. 

Recall that by Theorem 5.5, each pair of boundary curves of an Apollonian triangle 
have a common tangent line where they meet. Denoting the three such tangents 
the spline lines of the Apollonian triangle. Remark 4.4, and Lemma 4.3 give us the 
following: 

Lemma 6.1. The spline lines of an Apollonian triangle with vertices vi,V2,V3 are the 
median lines of the triangle AV1V2V3, and thus meet at a common point, which is the 
centroid of AV1V2V3 □. 

More crucially. Lemma 4.3 allows us to give a circle-free description of Apollonian 
curves. Indeed, letting c be the intersection point of the tangent lines to the endpoints 
Pi,P2 of an Apollonian curve 7, Lemma 4.3 implies (via Definition 5.1 and Theorem 
5.5) that the splitting point s of 7 is the intersection of the medians from pi,P2 of the 
triangle ApiP2C, and thus the centroid of the triangle ApiP2C. The tangent line to 7 at 
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s is parallel thus, by Observations 5.2 and 5.3, the following recursive procedure 
determines a dense set of points on the curve 7 given the triple (pi ,p2,c): 

(1) find the splitting point s as the centroid of Apip2C. 

(2) compute the intersections ci,C2 of the pic and p2C, respectively, with the line 
through s parallel to piP2- 

(3) carry out this procedure on the triples (pi,ci, s) and (5,02,^2)- 

By recalling that the centroid of a triangle lies 2/3 of the way along each median, 
the correctness of this procedure thus implies that the "generalized quadratic Bezier 
curves" with constant | described by Paoletti in his thesis [17] are Apollonian curves. 
Combined with Lemma 6.1, this procedure also gives a way of enumerating barycentric 
coordinates for a dense set of points on each of the boundary curves of an Apollonian 
triangle, in terms of its 3 vertices. Thus, in particular, all Apollonian triangles are 
equivalent under affine transformations. Conversely, since Proposition 4.1 implies that 
any 3 vertices in general position have a corresponding Apollonian triangle, the affine 
image of any Apollonian triangle must also be an Apollonian triangle. In particular: 

Theorem 6.2. For any three vertices vi,V2,V3 in general position, there is a unique 
Apollonian triangle whose vertices are vi,V2tV3. □ 

Another consequence of the affine equivalence of Apollonian triangles is conformal 
inequivalence of Apollonian triangulations: suppose 95 : 5 — > 5' is a conformal map 
between Apollonian triangulations which preserves the incidence structure. Let T and 
T' be their central Apollonian triangles, and a : T T' the corresponding affine map. 
By Remark 5.9, the points on dT computed by the recursive procedure above are the 
points at which T is incident to other Apollonian triangles of S; thus, cp — a on a, dense 
subset of dT, and therefore on all of dT- Since the real and imaginary parts of ip and a 
are harmonic, the maximum principle implies that if = a on T, and therefore on S as 
well, giving that S and S' are equivalent under a Euclidean similarity transformation. 
We stress that in general, even though T and T' are affinely equivalent, nonsimilar 
triangulations are not affinely equivalent, as can be easily be verified by hand. 

It is now easy to see from the right-angle requirement for Apollonian triangulations 
that the Apollonian triangulation associated to a particular triple of circles must be 
also be unique up to translation and scaling: by Remark 5.9, the initial choice of 
translation and scaling of the three degenerate Apollonian triangles determines the rest 
of the figure. (On the other hand, it is not at all obvious that Apollonian triangulations 
exist. This is proved in Theorem 7.1 below.) Hence by Proposition 4.1, an Apollonian 
triangulation is uniquely determined by the three pairwise intersection points of its 
three degenerate triangles: 

Theorem 6.3. For any three vertices vi,V2,v^, there is at most one Apollonian tri- 
angulation for which the set of vertices of its three degenerate Apollonian triangles is 

{vi,V2,V3}. □ 

To ensure that our piecewise-quadratic constructions are well-defined on a convex 
set, we will need to know something about the area of Apollonian triangles. Affine 
equivalence implies that there is a constant C such that the area of any Apollonian 
triangle is equal to C ■ A{T) where T is the Euclidean triangle with the same 3 vertices. 
In fact we can determine this constant exactly: 

Lemma 6.4. An Apollonian triangle T with vertices pi,p2,P3 has area ^A{T) where 
A[T) is the area of the triangle T — ApiP2P3- 
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Proof. Lemma 6.1 implies that the sphne hnes of T meet at the centroid c of T. It 
suffices to show that A{T Ci ApiPjc) — ^A{ApiPjc) for each C {1,2,3}; thus, 

without loss of generality, we will show that this holds ioi i = 1, j = 2. 

Let Ts = T n Apip2C, and let = Ap\pic \ Ts. We aim to compute the area of 
the complement using our recursive description of Apollonian curves. Step 1 of 
each stage of the recursive description computes a splitting point s' relative to points 
p'i,p'2tC' , and is the union of the triangles Ap'^p^s' for all such triples of points 
encountered in the procedure. As the median lines of any triangle divide it into 6 
regions of equal area, we have for each such triple that A{p'ip'2s') ~ ^A{p\p'2c') . 

Meanwhile, step 2 each each stage of the recursive construction computes new in- 
tersection points c'l , c'2 with which to carry out the procedure recursively. The sum of 
the area of the two triangles Ap'j^, c'j^, s' and As'c2P2 is 

A(Apl,c'i,5') + A(As',4,p'2) - ^-A{Ap\p'2s') - ^-A{Ap\p'2s') = ^A(Aplp^s'), 

Since ^A{Ap\p'2s') is the portion of the area of the triangle p'ip'2s' which lies between 
the lines piP2 and c'^, C2. Thus, the area A{Tf) is given by 

^(A.,«c) ■ (5 + (5) S + (5)1 + (5)1 + ^ ^ ■ ) ^ ?-^<^'"«* ° 

We conclude this section with some geometric bounds on Apollonian triangles. The 
following Observation is easily deduced from part 4 of Proposition 4.5: 

Observation 6.5. Given a proper Apollonian triangle with vertices Vi,V2, W3 generated 
from a non-initial circle C and parent triple of circles (Ci, C2, C3), the angles AviVjVk 
{{i,j,k} = {1,2,3}) are all > arctan(3/4) > ^ if C has smaller radius than each of 

Ci,C2,C3. □ 

Recall that Theorem 5.5 implies that pairs of boundary curves of an Apollonian 
triangle meet their common vertex at a common angle, and that there is thus a unique 
line tangent to both curves through their common vertex. We call such lines Li, L2, L3 
for each vertex vi,V2, fa the median lines of the Apollonian triangle, motivated by the 
fact that Lemma 4.3 implies that they are median lines of the triangle At;ii;2W3. 

Observation 6.6. The pairwise interior angles of the median lines Li,L2,-^3 of a 
proper Apollonian triangle all lie in the interval (^, ^). 

Proof. Part 1 of Proposition 4.5 gives that the interior angles of the median lines of 
the corresponding Apollonian triangle must satisfy < tt — |^ = |7r. The lower bound 
follows from ai + a2 + = 27r. □ 

7. Fractal solutions to the Apollonian PDE 

Our goal now is to prove that Apollonian triangulations exist, and that they support 
piecewise quadratic solutions to the Apollonian PDE which have constant Hessian on 
each Apollonian triangle. We prove the following theorems in this section: 

Theorem 7.1. To any mutually externally tangent circles Ci, C2, C3 in an Apollonian 
circle packing A, there exists a corresponding Apollonian triangulation S. Moreover, 
the closure of S is convex. 
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Theorem 7.2. For any Apollonian triangulation S there is a piecewise quadratic C^'^ 
map u : S ^ R such that for each Apollonian triangle Tc comprising S, the Hessian 
D^u is constant and equal to m{C) in the interior ofTc- 

Theorem 7.2 implies Theorem 1.3 from the Introduction via Proposition 2.8, by 
taking Y = S and Z = S, where S = 5(^1,^12,^3) is the Apollonian triangulation 
generated by the triple of circles c{Ai) for i = 1,2,3. Using the fact that S has 
full measure in S, proved in Section 7.2, this theorem constructs piecewise-quadratic 
solutions to the Apollonian PDE via Proposition 2.8. 

We will prove Theorems 7.1 and 7.2 in tandem; perhaps surprisingly, we do not see a 
simple geometric proof of Theorem 7.1, and instead, in the course of proving Theorem 
7.2, will prove that certain piecewise-quadratic approximations to u exist and use 
constraints on such constructions to achieve a recursive construction of approximations 
to S. 

7.1. The recursive construction. We begin our construction of u — and, simulta- 
neously S, which will be the limit set of the support of the approximations to u we 
construct — by considering the three initial matrices Ai = m{Ci) for i = 1, 2,3. 
Observation 2.6 implies that there are vectors vi^V2,v^ such that 

Ai = A^ + Vi® Vi for each i — 1,2, 3. 

We may then select distinct pi,p2,pz G such that Vi ■ (pj — Pk) = for {i,j,k} = 
{1,2,3}. Observation 2.5 and Definition 5.1 imply that we can choose degenerate 
Apollonian triangles 7^; corresponding to {Ai,Aj,Ak) {{i,j,k} = {1,2,3}) meeting 
at the points pi,P2,P3- Note that the straight sides of distinct Ta; meet only at right 
angles. 

It is easy to build a piecewise quadratic map uq e C^'^{Tai U 7a2 U Tas) whose 
Hessian lies in the set {^1,^2,^3}: for example, we can simply define uq as 

uo{x) := ^x*A^x + ^{v^ ■ {x - pj))'^ for x € Ta^ and i 7^ j. (7.1) 

We now extend this map to the full Apollonian triangulation by recursively choosing 
quadratic maps on successor Apollonian triangles that are compatible with the previous 
pieces. The result is a piecewise-quadratic C^'^ map whose pieces form a full measure 
subset of a compact set. By a quadratic function on we will mean a function of the 
form (p{x) = x*Ax + b* -x-l-c for some matrix A € 5*2, vector 6 e and c G M. Letting 
(1, 2, 3)^ denote {(1, 2, 3), (2, 3, 1), (3, 1, 2)}, the heart of the recursion is the following 
claim, illustrated in Figure 5. 

Claim. Suppose Bq is the successor of a triple {Bi,B2,B^), and that for each {i, j, k) e 
(1,2,3)^, we have that 7j is an Apollonian curve for (B^, Bj, B^) from p^ to pj, ipi is 
a quadratic function with Hessian Bi , and the value and gradient of ipi , (pj agree at pk 
for each k. 

Then there is a quadratic function ipo with Hessian Bq whose value and gradient 
agree with that of ipi at each qi :— 7i(ai(i?o)), and for each {i,j, k) e (1, 2, 3)^^, there is 
an Apollonian curve 7^ from qj to qk corresponding to the triple (Bq, Bj, Bk). (Here, 
the ai denotes the angle function a defined with respect to Bi.) 

We will first see how the claim allows the construction to work. Defining the level 
of each Ai,A2, A3 to be £{Ai) = 0, and recursively setting the level of a successor of a 
triple {Ai, Aj, Ak) as ma.x{£{Ai),£{Aj),£{Ak)) + 1, allows us to define a level-k partial 
Apollonian triangulation which will be the domain of our iterative constructions. 
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Pi 



Figure 5. The Apollonian curves 7^,7^' {i = 1,2,3) in the claim. 

Definition 7.3. A level-k partial Apollonian triangulation corresponding to {^1, A2, A3} 
is the subset Sk C 5(^1,^2,^3) consisting of the union of the Apollonian triangles 
Ta e 5 for which £{A) < k. 

Note that Uq is defined on a level-0 partial Apollonian triangulation. 

Consider now a C^'^ piecewise-quadratic function Ufc^i defined on the union of 
a level-(fc — 1) partial Apollonian triangulation Sk~i, whose Hessian on each G 
Sk-i is the matrix Ai. Any three pairwise intersecting triangles TAi,TAj ^Ta^ G ^fc-i 
bound some region R, and, denoting by 7^ the boundary curve of each Ta^ which 
coincides with the boundary of R and by Ps the shared endpoint of ^t^lu {{s,t,u} = 
{z,j,fc}), the hypotheses of the Claim are satisfied for {Bi, B2, B3) = {Ai,Aj,Ak), 
where ipi,ip2,(p3 are the quadratic extensions to the whole plane of the restrictions 
Uk-i\TA^ : Uk~i\rAj ' Uk-i\rA^ : respectively. 

Noting that the three Apollonian curves given by the claim bound an Apollonian 
triangle corresponding to the triple {Ai,Aj,Ak), the claim allows us to extend Uk-i 
to a C^'^ function Uk on the level-fc partial fractal Sk by setting Uk — ifio on the 
triangle Tac S Sk for the successor Ag of {Ai,Aj,Ak), for each externally tangent 
triple {Ai,Aj,Ak} in Sk-i- Letting U denote the topological closure of S, we can 
extend the limit u : 5 — >■ K of the Uk to a C^'^ function m : J7 — >■ K; to prove Theorems 
7.1 and 7.2, it remains to prove the Claim, and that 5 is a full-measure subset of its 
convex closure Z, so that in fact U = Z. We will prove that S is full-measure in 
Z in Section 7.2, and so turn our attention to proving the Claim. We make use of 
the following two technical lemmas for this purpose, whose proofs we postpone until 
Section 7.3. 

Lemma 7.4. Letpi,p2,P3 G be in general position, and let Vi — [pj — Pk)^ be the 
perpendicular vector for which the ray pi + svi (s G M+) intersects the segment pjpk, 
for each {i,j, k} — {1, 2, 3}. // ipi, ip2, y^a are quadratic functions satisfying 

D'^ipi = Ai, Dipi{pk) = Dipj{pk), and (pi{pk) = ^j{Pk) 

for each k} = {1,2,3}, where Ai = B~ + Vi Vi and tr(wi Vi) > 2(tr(_B) — 2) 
for some matrix B and vectors Vi perpendicular to pj —pk for each {i,j, k} — {1, 2, 3}, 
then there is a (unique) choice of Xq G K^, = Xq + t^Vi for ti/{vi ■ pj) > 1, and 
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6 e M^, c e M such that the map 




for k} — {1, 2, 3}. There is a matrix B and coefficients G M such that 



fori = 1,2,3. 

Observe now that in the settmg of the claim, the conditions of Lemma 7.4 are 
satisfied for Ai := Bi [i = 1,2,3), B := Bq and where Vi is the vector for which 
Bi — Bq = Vi ^ Vi for each i = 1,2,3; indeed Observations 2.5 and the definition of 
Apollonian curve ensure that is perpendicular to Pj —pk for each {i, j, fc} = {1, 2, 3}. 
Let now Xq, ti, and yi be as given in the Lemma. We wish to show that yi = "/i{a{Bo)) 
for each i. Letting Bij denote the successor of (Bq, Bi, Bj) for {i,j} C {1,2,3}, we 
apply Lemma 7.5 to the triples {pi,pj,yk} of points and {(pi,ipj,ipo} of functions for 
each of the three pairs {i,j} C {1, 2, 3}. In each case, we are given some matrix B for 
which 



for real numbers ak^i G K. 

Observation 2.6 now implies that either B = Bij or B = Bk; the latter possibility 
cannot happen, however: if we had B = Bk, then as p{Bk — Bq) = —p{Bo — Bk), 
Observation 2.5 would imply that yi ^yj is perpendicular to pi —pj. This is impossible 
since the constraint ts/(vs ■ pt) > 1 for {s,t} = {i,j} in Lemma 7.4 implies that the 
segment yiyj must intersect the segments PiPk and pjPk, yet part 1 of Proposition 4.5 
implies that ApiPjPk is acute. So we have indeed that the matrix B given by the 
apphcations of Lemma 7.5 to the triple {Bo,Bi,Bj) is Bij, for each {i,j} C {1,2,3}. 

For each {i,j,k} = {1,2,3}, Observation 2.5, Definition 5.1, Theorem 5.4, and the 
constraints (7.3), (7.4) now imply that yi = qi := ji{ai{Bo)), as the point 7i(ai(i3o)) 
is determined by the endpoints ^i{ak{Bj)),^i{ai{Bk)) and the condition from Def- 
inition 5.1 that ^i{ai{Bj)) - "fkiaiiBo)) and ji{ai{Bk)) - ji{ai{Bo)) are multiples 

1 /2 1 /2 

of Vi' {si{Bj,Bo)) and {si{Bk,Bo)), respectively (and so of pk - yi and pk - yj, 
respectively, by (7.3) and (7.4)). 

1 /2 

Similarly, the constraint (7.5) implies that qi — qj is a multiple of Vq (Bij) for the 

1/2 

function Vq defined with respect to the circle Bq. Definition 5.1 and Theorem 5.4 
now imply the existence of the curve 7^, completing the proof of the claim. 



D'^Lp, = B + a,(pj - pk)^ ® {pj - Pk) 



(7.2) 



B, = B^+ ak,s{Pk - yi) ® [Pk - yi) , 

Bj = B^ + ak^siPk - Vj)^ ® {Pk - Vj)^, and 

Bq = B^ + afe,o(y» - %)^ «) {yi - 2/j)^ 



(7.3) 
(7.4) 
(7.5) 



7.2. Full measure. We begin by noting a simple fact about triangle geometry, easily 
deduced by applying a similarity transformation to the fixed case of L = 1: 
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a 



Figure 6. To show that S has full measure in Z, we show that each 
Apollonian triangle Vi has area which is a universal positive constant 
fraction of the area of the region Ri it subdivides. Here, the boundaries 
of i?£ and are shown in long- and short-dashed lines, respectively. 



Proposition 7.6. Any angle a determines constants Ca,Da such that any triangle A 
which has an angle 9 > a and opposite side length £ < L has area A{A) < CaL'^, and 
any triangle which has angles Oi > a,92 > a sharing a side of length £ > L has area 

>DaL\ □ 

We wish to show that the interior of S has full measure in Z, defined as the convex 
closure of S. Recall that the straight sides of each pair of incident degenerate Apollo- 
nian triangles Vi,Vj {{i,j} C {1,2,3} intersect at right angles, so the 6 straight sides 
of Vi, V2, V^3 will form a convex boundary for Z. 

Letting thus = B \ St-, we have that is a disjoint union of some open sets 
Ri bordered by three pairwise intersecting Apollonian triangles, and Xt+i contains in 
each such region an Apollonian triangle Ve dividing the region further. To prove that 
the interior of S has full measure in Z, it thus suffices to show that the area A{V£) is 
at least a universal positive constant fraction k of the area A{Re) for each £, giving 
then that n(Yt) < (1 - ny-^niYi) -^t 0. 

For Re bordered by Apollonian triangles Vi,Vj,Vk and letting A' — AxiXjXk be 
the triangle whose vertices Xs are the points of pairwise intersections Vt , Ki for each 
{s,t,u} = {i,j,k} of the Apollonian triangles bordering Re, we will begin by noting 
that there is an absolute positive constant k' such that ^(A') > K'fj,{Re). For each 
{s,t,u} = {i,j,k}, the segment XgXt together with the lines L" and tangent to 
the boundary of Vu at Xs and Xt, respectively, form a triangle A^ such that Re C 
A' U Ai U Aj U Afe. Observations 6.5, 6.6, and 7.6 now imply that area of each A^ 
is universally bounded relative to the area of A', giving the existence k' satisfying 
A(A') > K'^i{Re). 

It thus remains to show that the Apollonian triangle Ve which subdivides £ satisfies 
f^iVe) > k"A{A') for some k" . (It can in fact be shown that /x(V^) = ^A(A') exactly. 
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but a lower bound suffices for our purposes.) Considering the triangle A" — Aabc 
whose vertices are the three vertices of Vi, there are three triangular components of 
A' lying outside of A"; denote them by A^, A^ , A'^ where A^ includes the vertex Xg 
for each s = k. The bound ZxgXtXu > f for each {s, t, u} = {«, j, k} together with 
Observation 6.5 implies there is a universal constant bounding the ratio of the area of 
A'g to A" for each s = k. Thus we have that the area of A' is universally bounded 
by a positive constant fraction of the area of A", and thus via Lemma 6.4 we have 
that there is a universal constant k" such that fi{Ve) > k"A{A'). 

Taking k — k' ■ k" we have that ^(Vf) > k^{Ri) for all £, as desired, giving that the 
measures ^J.{Yt) satisfy ^liYt) = (1 — k)*^^^(Yi) — >■ 0, so that ^{S) = n{Z). 

7.3. Proofs of two Lemmas. 

Proof of Lemma 7.4- We have 

(fi{x) — AiX + di — [B^ + Vi® Vi)x + di 

and the agreement of DLpj^ Dipk at Xi ({i, j, k} — {1, 2, 3}) together with Vi ■ {xj —Xk) — 
implies that the Dipi — B^x is constant independent of i on the triangle lSxiX2Xz: 
giving that 

Dip.i{x) = AiX — [vi® Vi)xj + d, (7.6) 

for a constant d G independent of i. Similarly, the value agreement constraints give 
that 

ipi{x) = -x^AiX — -x*j{vi (E) Vi)xj + dx + c, (7.7) 

for a constant c independent of i. Thus, by setting D := d adjusting C in the definition 
of Ml as necessary, we may assume that in fact c and d are 0. 

Fixing any point Xq inside the triangle X1X2X3, we define for each i a ray Ri ema- 
nating from Xq coincident with the line {tvi : t E Our goal is now to choose Xq 
such that there are points yi on each of the rays Ri satisfying the constraints of the 
Lemma. 

On each ray Ri, we can parametrize (pi := (pi{x) — ^x*B~x as functions fi{ti) = ^a^i^ 
{i = 1, 2, 3), and the function ifo := (po{x) — \x*'B~x as gi{ti) = ^b{t + hi)"^ + C , where 
ti is the distance from the line XjXk to x e Ri, hi is the distance from Xq to the line 
XjXk, and a,j and b are tT{viiS)Vi) and tr(_B), respectively. Moreover, since the gradients 
of (pi and ipo can both be expressed as multiple of Vi along the whole ray Ri, we have 
for any point x on Ri D Ui at distance t from XjXk that fi{ti) — g[{ti) implies that 
Dipi(x) = Dipo(x). Thus to prove the Lemma, it suffices to show that there are Xq 
and C such that for the resulting values of hi, the systems 

MU) = 9^{U) ^^^^ . r \a,tl = \b{U + h^f + C 

mu) = gm ™ explicitly, I ^ ^^^^ ^ 

have a solution over the real numbers for each i. 

It is now easy to solve these systems in terms of C; for each i, 

bh, 

ti = and hi - 



1 b-^ 

2 \" a,-b 



gives the unique solution. Note that > 2b ensures that the denominator in the 
expressions for hi (and ti) is positive for each i. Since \/—C takes on all positive real 
numbers and ti (Ai) — ai ~ b, there is a (negative) value C for which the distances hi 
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are the distances from the Unes XjXk to a point Xq inside AX1X2X3] it is the point 



/ tr(Ai)-tr(_B) \ ' 



With trihnear coordinates 

l<i<3 

The Lemma is now satisfied for this choice of C and Xq and for the points i/i on Ri 
at distance + ti from X^ for i = 1, 2, 3. □ 

Proof of Lemma 7.5. Let qi — p^—p2, 92 ~ Pi^Pa, and (73 = P2—P1 and := D^ipi. 
Since for any individual i — 1,2,3 we could assume without loss of generality that 
ipi = 0, the compatability conditions with ipj, ipk give 

UA,-A,)qj + iAk-A,)qk^O, 
[qjiAj - A,)qj - qiiAk - A,)qk - 

in each case. If we left multiply the first by ql. and add it to the second, we obtain 

ql{A,-A,)q,^0. (7.9) 

Since qi ■ qj 7^ 0, there are unique aij, f3ij, e M such that 

Aj -A,^ Uj^qf ® qf + pj^qf- ® q^ + 7,^(9,^ ® qf + qf (g) g,^), 

where {x,y)-^ = {—y,x). Since symmetry implies /3ji — —aij and (7.9) implies = 0, 
we in fact have 

A J - Ai = aj^qf ® qf ~ a^^qf «) qf. 
If we substitute this into (7.8), we obtain 

-a^jiqf ■ qj)q^ - a,k{qi ■ qk)qi = 0. 
Since qf ■ qj ~ —qf ■ qk, we obtain aij — aik. Thus there are G M such that 

At - A J = a^qf (g) qf ~ ajqf qf. 

In particular, we see that Ai — aiqf (g) qf is constant. □ 

7.4. Proof of Corollary 1.4. This is now an easy consequence of Theorem 7.2 and 
the viscosity theory, via Proposition 2.7. 

Proof of Corollary 1.4- Write v = Voc- Continuity of the derivative and value of v in 
U1UU2U U3 imply that 

v\u^ = ^x^A^x + Dx + C for i = 1,2,3 

for some D e M^, C G M. Let /3i be the portion of the boundary of R between Xj and 
Xk which does not include Xi, and let Vi be the vector perpendicular to Xj — x^ such 
that Xi + tvi intersects the segment XjX^. 

We let Vi = (it + tVi for t > 0. The V^^'s are pairwise disjoint. Thus, by first 
restricting the quadratic pieces J7i, C/2, C^3 of the map v to their intersection with the 
respective sets 1^, and then extending the quadratic pieces to the full V^'s, we may 
assume that Ui — Vi for each i = 1,2,3. 

We apply Lemma 7.5 to w|vi, wivi' ^Ivs! by Observation 2.6 there are up to two 
possibilities for the matrix B from (7.2); the fact that Ax-^X2Xj, is a acute, however, 
implies that we have that B is the successor of ^1,^2,^3- Thus letting S denote the 
Apollonian triangulation determined by a;i,a;2,a;3. Theorem 7.2 ensures the existence 
of a C^'^ map u which is piecewise quadratic whose quadratic pieces have domains 
forming the Apollonian triangulation S determined by the vertices Xi,a:2,X3. Letting 
U'i denote the degenerate Apollonian triangle in S intersecting Xj and Xk for each 
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Figure 7. The graph of the function c g C(IR^) over the rectangle 
[0,6] X [0,6/-\/3], where c describes the boundary 

dVtri = {lM{a,b,c{a,b)) -.a.beR} 

of Ftri- White and black correspond to c = 3 and c — A, respectively. 

{i,j,k} = {1,2,3}, we can extend u to a map u by extending the three degenerate 
pieces C// of S to sets V- = {x + tvi : x £ U[, t > 0}. Now we can find curves 7i from 
Xj to Xk lying inside VtHV-, and, letting f2 be the open region bounded by the curves 
7i; 72, 73, Proposition 2.7 implies that u + Dx + C and v are equal in fi, as they agree 
on the boundary = 7i U 72 U 73 . □ 

8. Further Questions 

Our results suggest a number of interesting questions beyond Conjecture 1.2. To 
highlight just a few, one direction comes from the natural extension of both the sandpile 
dymanics and the definition of F to other lattices. 

Problem 1. For the triangular lattice £tri Q generated by (1,0) and (1/2, -\/3/2), 
we define F-(-j.j to be the set of 2 x 2 real symmetric matrices A such that there exists 
u : Ctrl Z satisfying 

u > ^x^Ax and A^u < 5. (8.1) 

The algorithm from Section 3 can still be used in this case and we display its output in 
Figure 7. While the Apollonian structure of the rectangular case is missing, there does 
seem to be a set Vtri of isolated "peaks" such that ftri — What is the structure 
of these peaks? What about other lattices or graphs? To what extent does the set F 
characterize the underlying lattice/graph? 

Although we have explored several aspects of the geometry of Apollonian triangu- 
lations, many natural questions remain. For example: 

Problem 2. Is there a closed-form characterization of Apollonian curves? 
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Apollonian triangulations thcmsleves present some obvious targets, such as the deter- 
mination of their HausdorfF dimension. 
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